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Abstract. We develop an exact determinantal formula for the probability that the Airy2 
process is bounded by a function g on a finite interval. As an application, we provide a 
. . . direct proof that sup(.42(a;) — x 2 ) is distributed as a GOE random variable. Both the 

continuum formula and the GOE result have applications in the study of the end point of 
' _ _| , an unconstrained directed polymer in a disordered environment. We explain Johansson's 

■ [Joh03] observation that the GOE result follows from this polymer interpretation and 

exact results within that field. In a companion paper [MQR11] these continuum statistics 
. are used to compute the distribution of the endpoint of directed polymers. 

fC) ■ 1. Introduction 



Oh 
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The Airy2 process A 2 was introduced in Prahofer and Spohn [PS02] in the study of 
the scaling limit of a discrete polynuclear growth (PNG) model. It is expected to govern 
the asymptotic spatial fluctuations in a wide variety of random growth models on a one 
dimensional substrate with curved initial conditions, and the point-to-point free energies 
of directed random polymers in 1 + 1 dimensions (the KPZ universality class). It also 
arises as the scaling limit of the top eigenvalue in Dyson's Brownian motion [Dys62] for 
the Gaussian Unitary Ensemble (GUE) of random matrix theory (see [AGZ10] for more 
details). 

A 2 is defined through its finite-dimensional distributions, which are given by a Fredholm 
determinant formula: given xo, ■ ■ ■ , x n € R and to < • • • < t n in R, 

(1.1) PCA 2 (to) <x ,...,A 2 (t n ) <x n ) = det(/-f 1 / 2 K cxt f 1 / 2 ) L2({i0i ... A}xR) , 

where we have counting measure on {to, • • • , t n } an d Lebesgue measure on R, f is defined 
on {to, • • • , t n } x R by i(tj,x) = l x e(xj,oo)i an d the extended Airy kernel [PS02; FNH99; 
Mac94] is defined by 

K a ,, / ^_// OO ^e^')Ai(e + A)Ai(e' + A), if t > t> 
X : ^ ' U ~ [j^ d\ e-^') Ai(£ + A) Ai(£' + A), if t < f, 

where Ai(-) is the Airy function. In particular, the one point distribution of A2 is given by 
the Tracy- Widom largest eigenvalue distribution for GUE. 

K. Johansson [Joh03] proved the remarkable fact that 
Theorem 1. For every m G R, 



Pf sup (A 2 (t) - t 2 ) < m J = F GOE (4 1/3 m) 
VieR / 



Here -Fgoe denotes the Tracy- Widom largest eigenvalue distribution for the Gaussian 
Orthogonal Ensemble (GOE) [TW96]. It also arises as the one point distribution of the 
Airyi process, which governs the asymptotic spatial fluctuations in one dimensional random 
growth models with flat initial conditions, and the point-to- line free energies of directed 
random polymers in 1 + 1 dimensions. 



^The factor 4 1//3 corrects a minor mistake in Johansson's statement. See Section 2 for a discussion. 
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The proof of Theorem 1 in [Joh03] is indirect, using a functional limit theorem for 
the convergence of the PNG model to the Airy2 process, together with the connection 
between the PNG process and a certain last passage percolation model for which Baik and 
Rains [BR01] had proved the connection with GOE. In this article we develop a method 
to compute continuum probabilities for the Airy2 process -which is to say, compute the 
probability that the sample paths of the Airy2 process lie below a given function on any 
finite interval. This is then used to provide a direct proof of Theorem 1 starting only from 
determinantal formulas. 

Theorem 1 reflects a universal behaviour seen in a large class of one dimensional systems 
(the KPZ universality class starting with flat initial conditions) and therefore has attracted 
quite a bit of interest at the physical level. Much of the recent work is on finite systems 
of N nonintersecting random walks, the so-called vicious walkers [Fis84]. [Fei09; RS10; 
RS11] obtain various expressions for the maximum and position of the maximum at the 
finite N level. [FMS11] uses non-rigorous methods from gauge theory to obtain the GOE 
distribution in the large N limit, and furthermore connect the problem to Yang-Mills 
theory. 



Our computation of continuum probabilities starts with the following (earlier) variant 
of (1.1) due to Prahofer and Spohn [PS02], 

(1.2) F(A2(to)<X ,...,A 2 {tn)<X n ) 



H is the Airy Hamiltonian H = —d 2 + x and P a denotes the projection onto the inter- 
val (—00, a]. Here, and in everything that follows, the determinant means the Fredholm 
determinant in the Hilbert space L 2 (R). The equivalence of (1.1) and (1.2) was derived 
formally in [PS02] and [PS11]. In fact there are some subtleties, because, for example, it 
is not apriori obvious that for s,t > 0, e~ sH can be applied to the image of P a e~ tH . See 
[QR12] for a discussion of the technical details. 

Remark 1.1. The shifted Airy functions are the generalized eigenfunctions of the Airy 
Hamiltonian, as HAi(x — A) = AAi(x — A). The Airy kernel Km is the projection of H 
onto its negative generalized eigenspace. This is seen by observing that if we define the 
operator A to be the Airy transform, Af(x) := J*^^ dz Ai(x — z)f(z), then Km = APqA*. 

Fix £ < r. Given g £ ff 1 ([£, r]) (i.e. both g and its derivative are in L 2 ([£, r])), define an 
operator ©™ r i acting on L 2 (K) as follows: 0™ r i/(") = u ( r , •)> where u(r, •) is the solution 
at time r of the boundary value problem 



The fact that this problem makes sense for g € H l {[£, r]) is easy to prove and can be seen 
from the proof of Proposition 3.2 below. By taking a fine mesh in t we obtain a continuum 
version of (1.2): 

Theorem 2. 



det(/ - K Ai + P Xo e^-^ H P Xl e^~^ H ■ ■ ■ P Xn e^^ H K A ^j 



where K^i is the Airy kernel 




dtu + Hu 
u(£, x) 
u(t, x) 



for x < g(t), t E (£,r 

f(x)l x <g(e) 

for x > g(t). 



(1.3) 
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An expression in terms of determinants of solution operators of boundary value problems 
may not seem very practical. But in fact one can give an explicit expression for the kernel 
of the operator , in terms of Brownian motion. Let b(s) denote a Brownian motion 
with diffusion coefficient 2. By the Feynman-Kac formula, 



u(r,x) = E m=x (f(b(r))e St b{s)ds l b(s) < g{s) on [t 
The linear potential is removed by a parabolic shift, 
9f M /(x) = E m=x (f(b(r))e-K b ^ ds l b(s) < g{s) on [M ) 

= E 6W=x (/(6(r))e^W-^«+^ 3 -^)/ 3 +/^ d ^)-/^ 2 ^l bW < 9(s) on ( , 
= E m=x _ P (f(b(r) + r 2 )e^W^)-( fe «^ 2 )+^-^)/ 3 l b(s)+s2 < 5(s) on M ) , 



where in the second equality we used integration by parts and added and subtracted (r 3 
£ 3 )/3 and in the third one we used the Cameron-Martin-Girsanov formula. This gives 

Theorem 3. Let Q?^ r Jx,y) denote the integral kernel o/Oj^,. Then 



\e)=x-e 2 Mr)= y -r 2 ^( s ) ^ 9(s) - s 2 on [£,r]J , 



where the probability is computed with respect to a Brownian bridge b(s) from x — £ 2 at 
time I to y — r 2 at time r and with diffusion coefficient 2. 

This gives a formula which can be used in applications. The obvious one is the case 
g(t) = t 2 + m, in which the probability can easily be computed by the reflection principle 
(method of images). A second one is the computation of the joint distribution of the max 
and argmax of the Airy2 process minus a parabola, which appears in a companion paper 
[MQR11]. The simple result in the case g(t) = t 2 + m, setting — I = r = L, is that 

where Rl is the reflection term 

(1.5) R L (X, y) = _^ e -{x+y-2™-2L*f/SL- {x+ y)L + 2L3/^ 

The first term in Ox, has been reexpressed in terms of the Airy Hamiltonian by reversing 
the use of the Cameron-Martin-Girsanov and Feynman-Kac formulas. 

To obtain the L — > 00 asymptotics, decompose ©x so as to expose the two limiting 
terms, as well remainder term fir • 

(1.6) L = e~ 2LH - R L + n L , 

where Ox = (Rl — P m +L 2 RLPm+L 2 ) ~ (e~ 2LH — P m+ L 2e ~ 2LH Pm+L 2 ) ■ I n Section 5 we will 
show that 

Lemma 1.2. As L goes to infinity, 



n L := e LH K Ai n L e LH K Ai -> 



in trace norm. 



Referring to (1.3), we have by the cyclic property of determinants and the identity 
e 2LH K Ai = (e LH K Ai ) 2 that 

(1.7) F(A 2 (t) < g(t) for t G [-L, L]) = det(l - K Ai + e LH K Ai @ L e LH K Ai ) . 



CONTINUUM STATISTICS OF THE AIRY 2 PROCESS 



1 



Since e Kpse e Kai = i^Ai an d due to Lemma 1.2, one sees that the key point is 
the limiting behaviour in L of e LH K aiRls LH K ai- Remarkably, it does not depend on L 
and gives the kernel of -Fgoe, thus providing a proof of Theorem 1. 

Proposition 1.3. For all L > 0, 

e LH K Ai R L e LH K Ai = AP RP A*, 
where the A is the Airy transform (see Remark 1.1), and 

R(X, A) := 2" 1 / 3 Ai(2~ 1 / 3 (2m - A - A)). 

Furthermore, 

(1.8) det(j - APqRPqA*} = F GOE (4 1/3 m). 

The last equality is a version of the determinantal formula for -F G oe proved by Ferrari 
and Spohn [FS05], and which essentially goes back to Sasamoto [Sas05]: 

F GO E(m) = det (J - P B m P ), where B m (x, y) = M(x + y + m). 

This can be seen as follows. Using the cyclic property of the determinant and the reflection 
operator o~f(x) = f(—x) we may rewrite the determinant in (1.8) as 

(1.9) det (7 - P RPq) = det (j - aP RP a^ = det (i - P aRaPo) , 

where we have used that AA* = a 2 = I. On the other hand we have aRa(X, A) = 
2~V3 Ai(2 _1 / 3 (A + A + 2m)). Performing the change of variables A i->- 2 1 / 3 A, A H> 2 1 / 3 A in 
the Fredholm determinant shows that the determinants in (1.9) equal det(/ — Po-B 4 i/3 m -Po)- 

The rest of the paper is organized as follows. In Section 2 we give an overview of the 
approach of Johansson [Joh03] explaining how Theorem 1 can be obtained indirectly using 
the connection of the Airy2 process with last passage percolation. Section 3 contains a 
brief introduction to relevant ideas of Fredholm determinants and then provides a proof 
of Theorem 2. Section 4 provides a short proof of Proposition 1.3. Finally, Section 5 
is devoted to the proof of Lemma 1.2, which essentially amounts to asymptotic analysis 
involving the Airy function. 

Acknowledgements. JQ and DR were supported by the Natural Science and Engineering 
Research Council of Canada, and DR was supported by a Fields-Ontario Postdoctoral Fel- 
lowship. IC was supported by NSF through the PIRE grant OISE-07-30136. The authors 
thank Victor Dotsenko and Konstantine Khanin for interesting and helpful discussions, and 
Kurt Johansson for several references to the physics literature. Part of this work was done 
during the Fields Institute program "Dynamics and Transport in Disordered Systems" and 
the authors would like to thank the Fields Institute for its hospitality. 

2. Indirect derivation of Theorem 1 through last passage percolation 

As we mentioned in the introduction, Johansson [Joh03] presented an indirect proof of 
Theorem 1 by way of the PNG model. His idea was entirely correct, but in the process of 
translating between the available results at the time, a factor of 4 1//3 was lost. The purpose 
of this section is to explain Johansson's approach and account for the missing 4 1 / 3 . 

We consider the PNG model (which we define below) with two types of initial conditions 
(droplet and flat), and show that by coupling them to the same Poisson point process 
environment we can represent the one-point distribution for the flat case as the maximum 
of the interface in the droplet case. Asymptotics of this relationship leads to the identity 
in Theorem 1. 
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(x, t) (x, t) (x, t) 




• g(x) = 



space space space 

Figure 1. The maximization problems coupled to the same Poissonian 
environment. Paths tt must be Lipschitz 1 functions of time and T(ir) 
represents a count of the number of Poisson points encountered by tt. Left: 
A general function g{x) represents the possible starting space-time starting 
location. Middle: The droplet geometry in which g(x) = \x\. Right: The 
flat geometry in which g(x) = 0. 



Consider a space-time Poisson point process P of intensity 2. Define a height function 
above x at time t as 

h g (x,t) = max T(ir) 

w.g^(x,t) 

where g represents a space-time curve (g(x), x) xG k, 7r is a Lipschitz 1 function of time (i.e., 
\tt(s) — tt(s')\ < \s — s'\ for all s,s'), ir : g — ?■ (x,t) means that tt starts at a point of the 
form (g(x),x) and ends at the point (x,t), and T(ir) represents the sum of the number 
of Poisson points that it touches. We will specialize this definition to two cases. In the 
droplet geometry (for which we write /j dro P let ) we take g = \x\, hence we only consider paths 
originating along a wedge. As a result the maximal path will always originate at the origin 
(0,0). In the flat geometry (for which we write /i flat ) we take g = 0, hence we consider 
Lipschitz paths starting in any spatial location at time and ending at x at time t. This 
are illustrated in Figure 1. 

Couple P to another Poisson point process P via P{A) = P(rtA), where for any Borel 
set A £ M 2 , (y, s) £ TtA if and only if (— y, t — s) £ A (one should think of this as a time- 
reversal of the Poisson point process where s i— > t — s and x h-> —x). Let /i flat represent 
the flat geometry height function built on the P Poisson point process. Then the following 
relation holds 

h^\t,Q) =max/i droplct (t,x). 

Asymptotic fluctuation statistics have been derived for both the droplet and flat geome- 
tries and (up to justification of taking the limit inside the maximum, as done in [Joh03] for a 
related model) the limiting statistics also respect the same relationship above. Specifically 
[PS02] (see also [BFP08] for the specific choices of scaling used below) shows that 

.. h d ™P lct (t,t 2 / 3 x)-2t 2 
lim = A2{x) - x . 

This implies that (up to the justifications mentioned above) 

,. /i flat (t,0) -2t , . . , 2 , 

lim — = max \A2\x) — x . 

t^oo £1/3 xGR ' 

On the other hand, [BFS08] shows that 

lim ^(*.0)-2 t = 

t->oo £l/3 1V ' 
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where A\ is the Airyi process. Combining these two identities shows that 

P(maxL4 2 (x) - x 2 ) < m) = PU(0) < 2~ 1 / 3 m) = F G o E (4 1/3 m), 

where the last equality follows from work of Ferrari and Spohn [FS05] which shows that 
P(.Ai(0) < m) = F GOE (2m). 

3. Proof of Theorem 2 

The operator in (1.2) should be seen as a discrete version of the boundary value problem 
operator 0^ In particular, for n > let t{ = £ + i(r — £)/n, i = 0, . . . , n, and define the 
discrete time boundary value problem operator 

The proof of Theorem 2 amounts to showing that, as n goes to infinity, the discrete operator 
converges to the limiting operator 0j^ , . This convergence must be in a suitably strong 
sense to ensure the convergence of the Fredholm determinants. Therefore, before turning 
to the proof of Theorem 2, let us briefly review some facts about Fredholm determinants, 
trace class operators and Hilbert-Schmidt operators (see Section 2.3 in [ACQ11] for more 
details, a complete treatment can be found in [Sim05]). 

Consider a separable Hilbert space T~L and let A be a bounded linear operator acting on 
T~L. Let \A\ = V A* A be the unique positive square root of the operator A* A. The trace 
norm of A is defined as ||^4||i = X^^i( e m l^l e n); where {e n } n >i is any orthonormal basis 
of %. We say that A £ Bi(T-L), the family of trace class operators, if ||A||i < oo. For 
A £ £>i(%), one can define the trace tv(A) = X^^i( e «> Ae n ) and then the Hilbert-Schmidt 
norm \\A\\2 = ^/tr(\A\ 2 ). We say that A £ B 2 ('H), the family Hilbert-Schmidt operators, if 
\\A\\2 < oo. The following lemma collects some results which we will need in the sequel, 
they can be found in Chapters 1-3 of [Sim05]: 

Lemma 3.1. 

(a) Ai4 det(/ + A) is a continuous function on B\{T-L). Explicitly, 

| det(J + A) - det(J + B)\ < \\A - B\\x exp(||A||i + + 1). 

(b) If A £ B X (H) and A = BC with B,C £ B 2 {H), then \\A\\ X < ||fi|| 2 ||C||2. 

(c) If ||j4||o P denotes the operator norm of A in H, then ||-A|| p < ||^4||2 < ll^lli? 
\\AB\\i < \\A\\ op \\B\\i and \\AB\\ 2 < ||A|| op ||-B|| 2 . 

(d) If A £ B 2 ('H) ) then \\A*\\ 2 = ||A|| 2 . If A has integral kernel A{x,y), then 

1/2 

\A\\ 2 = ( / dxdy\A{x,y)\ 2 



The proof of the continuum limit of (1.2) will follow easily from the next proposition. 
Proposition 3.2. Assume g £ H 1 ^, r]). Then the operators Kai — Q^Mr] H ^ Ai 
and K A \ - @f £ , r f {r ~ £)H K A i are in £i(L 2 (M)), with \\K M - Ql^^e^ 11 K A \\\i bounded 
uniformly in n. Furthermore, for any fixed i < r we have, writing = 2 k , 
(3.1) km ||(A Ai - ®l M e^) H K M ) - (K M - Q^./^ K M )\\i = 0- 

The idea of the proof is the following. Just as done in the introduction for 0^ , , 
it is possible to use the Feynman-Kac and Cameron-Martin-Girsanov formulas to write 
a formula for the kernel of Q 9 ,„ , in terms of a path integral with a killing potential 

enforced only at the dyadic mesh of times {tj}™^ (as opposed to being enforced at all 
times in [£, r] ) . If one considers a parabolic barrier g then the kernel for 0^ ^ , is given in 
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terms of the probability of a Brownian bridge exceeding a fixed value at some time {tj}™^. 
This is compared to the analogous kernel for Gj^ , given in terms of the probability of a 
Brownian bridge exceeding a fixed valued at any time t G [£, r]. As the mesh goes to zero, 
these two probabilities converge and hence so do the kernels. This proves the proposition 
for parabolic g, and the extension to g G r]) then follows readily since H 1 is the 

Cameron-Martin space for Brownian motion. 

Proof of Proposition 3. 2. We will first prove the result assuming that g(s) = {s-\{£+r) 2 ) 2 . 
Let (f(x) = (1 + x 2 ) 1 / 2 and define the multiplication operator Mf(x) = ip(x)f(x) (note 
that the choice of p is not particularly important and any strictly positive, polynomially 
growing function would do). To estimate the trace norm of K A \ — Gj^ r ^ r ~^ H Kai we use 
Lemma 3.1 to write 

(3.2) \\K M - Q^-^KmI < \\(e-^ H - Q^ r] )M\\ 2 \\M~ l e^ H K A] \\ r 
For the second Hilbert-Schmidt norm above we have by (4.1) that 



|| M -i e (r-^ KA .i|i / dxdy d\d\tp{x)- 2 e {x+x){r - e) Ai(x- A)Ai(y-A) 

Jr 2 J(-oo,0] 2 

(3.3) • Ai(x — A) Ai(y — A) 



00 



dx / dA^OrrV^^AiOr-A) 2 < c(2(r - ^))^ 1 ||^ 1 ||^, 



where c = max xe K Ai(x) 2 < 00. 

Now we consider the first norm on the right side of (3.2). Shifting time by —(£ + r)/2 
in the definition of @y r ] it is clear that this operator equals @®_ L L y where L = (r — £)/2 
and g(s) = s 2 . Using the formula for @^_ L n(^> y) given in Theorem 3 we get 

e -(x-y) 2 /8L-{x+y)L+2L 3 /3 , . 

%,rpiV) = ^ n ( -L)= X -L*,kL)=y-A b ( S ) < 011 H^l) ■ 

Similarly, the kernel of e~^ r ~^ H = e~ 2LH equals the above one with the probability replaced 
by 1, and hence 

-(x-y) 2 /8L-(x+y)L+2L 3 /3 

(3.4) (e-(r-l)H _ @9 [e r] ) M (x, y) = -j== V (y) 

L2 ( b(s) > for some s G [— L, L] ) . 



' r b(-L)=x-L 2 ,b(L)=y- 

Using a known Brownian bridge formula (see for example page 67 in [BS02]), the lat- 
ter crossing probability equals e~( x ~ L2 ^ y ~ L2 ^ 2L if x < L 2 ,y < L 2 and 1 otherwise, and 
therefore 

(3.5) 

Ue -(r-t)H _ Q9 )M\\ 2 = -*-[ dxdy [ e -(-^) 2 /8L-(, +y )L-4L3/3 ]2 ^ (2/ + ^2 

+ — I dxdy r e -(^) 2 /8L-(x +?/ )L-4L3/3l2^ (y + ^2 

8vrL yf_ 0O ,ni2 



where we have performed the change of variables x \— > x + L 2 , y 1— >• y + L 2 . Both Gaussian 
integrals can be easily seen to be finite, so we have shown that (e~( r ~^ H — Gj^ r j)M G 

B 2 (L 2 (R)). Using this with (3.2) and (3.3) it follows that K Ai - @ 9 [e r] e^ T '-^ H K Ai is in 

£i(L 2 (R)). 
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Next we observe that we can shift time and apply the Feynman-Kac and Cameron- 
Martin-Girsanov formulas directly on Q 9 n ^ (n times) exactly as we did for 6^ r p and it 
is not hard to check that we get a formula analogous to (3.4): 

( 6 -0-Ofl _ Q^ r] )M(x,y) = -±= e -^ ^-^ L+2L3 '^{y) 

■ IP- 

^ b n (-L)=x-L 2 ,b n (L)=y- 



L 2[b n {ti) > for some i G {0, . . . , n}\ . 



where b is now a discrete time random walk with Gaussian jumps with mean and variance 
2L/n, started at time — L at x — L 2 , conditioned to hit y — L 2 at time L, and jumping 
at times f" = — L + 2iL/n, i > 0. A simple coupling argument (see the next paragraph) 
shows that the last probability is less than the corresponding one for the Brownian bridge, 
and thus we obtain for ||i^Ai — u r ^ r ~^ H K^iWi the same bound as the one we get for 

II-Kai - Qf^f^^KAiWi from (3.5). This bound is, in particular, independent of n. 

Finally, in order to prove (3.1) we couple the Brownian bridge b and the conditioned 
random walk b Uk by simply letting b nk (t™ k ) = b(t™ k ) for each i = 0, ...,njfc. Since the 
Brownian bridge hits the positive half-line whenever the conditioned random walk does, it 
is clear that 



e -(x-y) 2 /SL-(x+y)L+2L 3 /3 



<-p(y)qn k (x,y), 



where q nk (x,y) is the probability that the Brownian bridge b(s) hits the positive half-line 
for s E [-L, L] but not for any s G {*Q fc , . . . , t^n }• Since every point is regular for one- 
dimensional Brownian motion, q nk (x,y) \ as k —> 00 for every fixed x,y, and thus by 
the monotone convergence theorem we deduce that \\{e~^ r ~^ H — Q^ r ^)M — {e.~^ r ~^ H — 
®n fc ^ r j)M||2 — > as k — > 00. Using (3.3) and a decomposition analogous to (3.2) yields 
(3J).' r 

To extend the result to g € H 1 ^, r]) we note that everything in the above argument deals 
with properties of a Brownian motion b(s) killed at the positive half-line. In the general 
case we will have by Theorem 3 a Brownian motion b(s) killed at the boundary g(s) — s 2 
or, equivalently, a process b(s) = b(s) — g(s) + s 2 killed at the positive half-line. Using 
the Cameron-Martin-Girsanov Theorem we can rewrite the probabilities for b(s) in terms 
of probabilities for b(s). Since g(s) is a deterministic function in H l ([l, r]), the Radon- 
Nikodym derivative of b(s) with respect to b(s) has finite second moment, and thus by using 
the Cauchy-Schwarz inequality we get the first two statements in the result from the above 
arguments. The convergence in (3.1) follows as well from the above arguments because it 
only depends on almost sure properties of the corresponding Brownian motion. □ 



Proof of Theorem 2. With the notation introduced before Proposition 3.2 we have 

PG4 2 (io) < gfo), . . .,A 2 (t n ) < g(t rik )) = det(j - K Ai + e 9 nk[er] e^ H K A ^j , 

where = 2 k . Since the Airy2 process has a version with continuous paths (see [Joh03; 
CH11; QR12]), the probability above converges to P(.A 2 (i) < g{t) for t G [£,r]) as k — > 00. 
The result now follows from Proposition 3.2 and Lemma 3.1, which imply that 

lim det(J - K M + Q 9 nM e^ H K M ) = det(J - K Ai + Q 9 [e r] e^ H K Ai ). □ 
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4. Proof of Proposition 1.3 

Since K A \ is the projection onto the negative (generalized) eigenspace of H (see Remark 
1.1), we have 

r° 

(4.1) e LH K M (x,z) = dX e XL Ai(x — A) Ai(z — A) dX. 

J — oo 

Then, recalling that the Airy transform is given by Af(x) = dXAi(x — X)f(X), we can 
write 

(4.2) e LH K Ai R L e LH K Ai = AP R L P A*, 
where 

£ L ( A , A) = —L= / dz dz e -(-+S-2— 2^)V8L-(,+5)L+(A+A)i+2^/3 Ai(z _ A) Ai(5 _ \). 
V8ttL Jr2 

Applying the change of variables 2u = z + z, 2v = z — z, we get 

R L (X, A) = -jL= [ dudv e - ( "- m 2 - L2)2 -2 M L+(A+A)L+§L3 ^ + y _ ^ ^ _ y _ ^ 
Using the formula 

/oo 
da; Ai(a + x) Ai(6 - x) = 2" 1/3 Ai(2" 1/3 (a + 6)) 
-oo 

(see, for example, (3.108) in [VS10]), the v integral equals 2" 1 / 3 Ai(2~ 1 / 3 (2n - A - A)). 
Therefore 

o— 1/3 roo 

Rl{X, A) = ^= / due -(-— i 2 ) 2 /2L-2«L + (A + A)L+2LV3 Ai(2 -l/3 (2u _ A - A)) 

y2nL J -oo 

o— 1/3 i /■ roo 

(it / ^ ne -(«-m-i 2 ) 2 /2i-2«L+(A+A)L+2L 3 /3+t :! /3-2- 1 / 3 t(A+A-2M) 



2vri J T 

J —oo 

where in the second equality we have used the contour integral representation of the Airy 
function, Ai(x) = i J r dt e* / 3 ~ tx with F = {c+is and c any positive real number. 

The u integral is just a Gaussian integral, and computing it we get 

R L (X,X) = 2 " 1/3 / ^ e t 3 /3+2V3Lt2 + [4V3 i 2 +2 -l/3 (A+ X_ 2m )]t + 2i3/3+(A+A)L_ 

2vri J r 

Now we perform the change of variables t = s — 2 1//3 L to obtain 

£ L (A,A) = — [ ^e s3 / 3 - 2 " 1/3 ( 2m - A -^ = 2- 1 /3Ai(2- 1 /3( 2m -A-A)) 

27T1 J T i 

(here the contour T' is simply T shifted by 2 1 / 3 L, so the integral still gives an Airy function). 
Note how all the terms involving L have canceled. 

5. Proof of Lemma 1.2 

The proof of this result amounts to asymptotic analysis of integrals involving the Airy 
function. The following well-known estimates for the Airy function (see (10.4.59-60) in 
[AS64]) go a long way in the proof: 

(5.1) |Ai(x)| < Ce"^ 372 for x > 0, |Ai(x)| < C for x < 0. 

However, at one important point the above bounds for x < will prove insufficient, and 
it will become necessary to utilize a representation (5.6) for Ai(x) which splits it into 
complex oscillations of opposite phase. Then, following standard methods of asymptotics 
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for oscillatory integrals (i.e., shifting real contours up and down to turn oscillations into 
exponential decay), we will achieve our bounds needed to complete the proof of this lemma. 

We will use the following version of Laplace's method, which we state without proof 
(see, for instance, [Erd56]): 

Lemma 5.1. Let 

I{M)= f dxf(x)e^ x)M , 



where C M. n is a (possibly unbounded) open polygonal domain and f and if are smooth 
functions defined on 0. Assume that the local maxima of ip are attained at a finite subset 
{xi, . . . , x n } of £1. Then there is a constant C > such that 

n 

\I{M)\ < Cj2 M ~ Ki \f( x i)\ e<p(Xi)M 

k=l 

for large enough M , where Ki = (n + l)/2 if xi £ d£l and Ki = n/2 if Xi £ Q. 

We write Cl^ = Q} L — where 

Vl\ = e LH K Ai (R L - P m+L 2R L P m+L 2) e LH K M , 

il L = e K A i{e - t> m+L 2e P m+L 2)e K Ai , 

The proof of Lemma 1.2 is contained in the next two lemmas. 
Lemma 5.2. 

"^LlU 



10} lli ► 0. 



L— >oo 

Proof. We proceed as in (4.2) and factorize £l l L as 

nl = ACl\A\ 

where 

(5.2) ni(X, A) = -jL= J_ dz d~Z e -(*+2-2rn-2L*)y8L-{z + S)L+2L*/3+(X + ~X)L 

■Ai(z- A)Ai(z- A)1 a X < , 

with D = R 2 \ (— 00, m + L 2 ] 2 . Using the Plancherel formula for the Airy transform 
J f 2 = f(Af) 2 , we have ||A|| op = ||yl*|| p = 1, so by Lemma 3.1 it will be enough to show 
that 

(5.3) ||f^||i >0. 

L— >oo 

Performing the change of variables z = L 2 (l — w) + m, z = L 2 (l — w) + m in (5.2) the 
kernel becomes 

(5.4) nl(\,\) = =— = / dw dw e f( w,w } Ai(L 2 (l — w) + m — \) 



V8tt Jd 

■Ai(L 2 (l-w)+m-~X)l xMo , 

where D = M. 2 \ [0, oo) 2 and f(w,w) = ~( w +^ 2 + ( w + w) - §. 

We split the region D into the union of three disjoint regions of pairs (w, w): D\ = {w < 
l,w < 1} - {0 < w < 1,0 < w < 1}, D 2 = {w < 0,w > 1} and D' 2 = {w > l,w < 0}. 
By the triangle inequality we can bound ||0^(A,A)||i by the sum of the trace norms of the 
operators obtained by restricting the integral in (5.4) to each of the regions D±, D 2 and 
D' 2 . We will write for the operator restricted to D\, with the analogous notation 
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for the other regions. Notice that, due to the symmetry of our formula, we do not need to 
bound ||n^l£)/||i, as it satisfies the same bound as ||n^l£) 2 ||i. 

Let us focus first on the operator restricted to D±, which is the simplest case because 
1 - w > and 1 - w > 0. We write Cl\ = GiG 2 G 3 with 

r 7/2 

G x {\,w) = e^ L Ai(L 2 (l -w)+m- X)l x < . w <i, G 2 (w, w) = —e L f ^' a h {w ^ )eDl , 

\/8ir 

G 3 (w, A) = e^ L Ai(L 2 (l -w)+m- ^)^x<o,w<i 

and use Lemma 3.1 to estimate || 1 < 1 1 Ci 1 1 2 1 1 C2 1 1 2 1 1 *-^3 1 1 2 • We begin with G\ and assume 
first that w < 1 + mi" 2 , so that using (5.1) have 

/0 /•min{l,l+"iL -2 } 
d\ I dw e 2{x ~ m ^ L Ai(L 2 (l - w) + m - A) 2 

-OO J —OO 

poo poo 

< Ce~ 2mL / dX / dwe- 2X -^ L2w+x ^ /2 < C'e- 2mL . 



Likewise, using the other bound in (5.1) it is easy to get ||GiP 1+m ^-2 1| 2 < C\m\e~ 2mL . 
Therefore ||Gi||2 < Ce~ mL , and of course the same bound works for G3. On the other 
hand, / attains its maximum on D\ at the points (1, 0) and (0, 1), where its value is — 14. 
Lemma 5.1 then allows to conclude that 



|G 2 || 2 = f dwdw—e 2L3 K w ^<CL 7 (L 3 r 3 / 2 e- u 
Jd 1 8tt 



L1L 3 /12 

Putting the three bounds together we deduce that 



(5.5) ll^ilDilli <e~ CL3 

for some C > 0. 

Let us now turn to the trace norm of the operator restricted to D2 (and hence also to D' 2 ). 
This bound is slightly harder owing to the fact that one of the Airy functions is oscillatory 
(rather than rapidly decaying) in this region. As readily derived from the contour integral 
representation of the Airy function by deforming the contour and performing a change of 
variables, Ai(-) may alternatively be expressed as 



Ai(x) = Re 



-x 



2vri 



/ ds exp(i(— x 



where T is the contour {s = a + b(a)i : a > 0} with b(a) = (a — g w^. This contour is 



3a 

,3, 



the steepest descent contour for f(s) = i(— s + s /3) and has the property that Im/(s) = 
Im/(so) = —2/3, where so = 1 is a critical point of /. Along T we can write f(s) = 
— 2/3i + g(s) where g(s) is real valued, g(so) = and g(s) decays to —00 monotonically 
and quadratically with respect to | s — so \ ■ Thus we may also write 



(5.6) Ai(x) = \{G{-x) + G{-x)) 

where 

G{x) =exp(-|x 3/2 i)^ f ds eMx 3/2 9(s)). 
2vri J r 

This expansion of the Airy function is the key to our oscillatory asymptotics. 

By applying the change of variables w = + L~ 3 / 2 t> and w = 4 + L~ 3 ^ 2 v the integral we 
wish to bound is given by 
(5.7) 

j 1/2 (A+A)L rO poo . .,2 

n\\,\) = = / dvh L (v) dve~^- AU-3L 2 - L 1/2 v + m - X)l, l<n , 

V8vr J-oo J-3L3/2 
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where 

We rewrite this as 

(5.8) 

with 



h L (v) 



fli(A,t;) 



e3 L ' 3 Ai(L 2 - L l ' 2 v 
Q} L = H1H2 

L l/2 



m — A). 



XL~v 2 /16 



dv e~ 



3L3/2 



- / 16 Mt>)1 A| «<o, 
1 Ai(-3L 2 - L 1/2 £ + m - X)l v ~ x 



We will focus on the last integral in v and prove that it is bounded by e~ CIj3 for some 
C > 0. As the Airy function is bounded on the real axis, we readily find that due to 
the Gaussian term, we may cut our integral outside of a region R$ = (—5L 3 ' 2 ,5L 3 ' 2 ) by 
introducing an error of order e 

XJi 
with 



-CL 3 ^ Thus we may restrict our attention to Rs- 
Using the expansion given by equation (5.6), the v integral may be written as + 



II 



dv e —G(3L 2 + L / v — m + X). 



5L :i / 2 



T7I ^ „-CL 3 



We wish to show that = < e 



For simplicity we set v = A 



m 



0, though 



the argument below does not rely on this assumption and applies equally well for all A < 0, 
v < and m / as necessary. Under this simplification, and performing a change of 
variables from v to r by setting 



L 3 r = (3L 2 + L 1 / 2 vf/ 2 , 



we obtain 



2 L 3/2 



3 3 / 2 +<5' 



3 3 / 2 -5' 



drr _ 1 /3 e -4( r 2/3_3 ) 2 G((L 3 r)2/3 ^ 



Since we can consider an arbitrary 5 before the change of variables, we can likewise consider 
an arbitrary 6' > for which to bound It,. Plugging in the expression for G we get 



II 



£5/2 



3 3 / 2 +<5' _ L 3 

dr e 

3 3 / 2 -<5' 



,2/3, 



exp(L 3 rg(s))ds. 



Observe that this integrand is analytic in r. Thus by Cauchy's theorem, rather than 
integrating from 3 3 / 2 - 6' to 3 3 / 2 + 5' along the real axis, we may do so along any other 
curve between these points. Due to the properties of g(s) along T, as long as Re(r) > we 
have that 



ds exp(L 3 r<?(s)) 



< 



ds exp(L 3 Ke(r)g(s)), 



which is certainly a bounded function of r for Re(r) > 0. The decay of the integrand is 
thus controlled by 



(5.9) 



Re 



Hr 2/3 ~ 3) : 



n 



Informed by this we may deform the r integration contour to the contour B = B1UB2UB3 
where B x = {3 3 / 2 - 6' + iy : y € [0,??]}, B 2 = {x + irj : x E [3 3 / 2 - <5',3 3 / 2 + <$']} and 
B3 = {3 3 / 2 + 5' + iy : y € [0, r/]}. It is an exercise in basic complex analysis to see that 
one can choose rj in such a way that, along the contour B, (5.9) stays bounded below a 
constant —C for C > 0. This implies that the exponential is bounded by e~ CL along that 
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curve and hence for the entire integral we get < e GL> for some C > 0. Going back to 
the definition of H 2 this implies that 



\H 2 \\ 2 = U dvd\e 2 ~ XL - v2 l\l L +T L f < e~ CL \ 
4 -'-(00, o] 2 



Returning to (5.8), it now suffices by Lemma 3.1 to prove that ||i?i||2 does not grow like 
e . This follows readily by integration using (5.1), which implies that 

(5.10) lGg(M«)) ~ - l(L 2 - L^v + Xf 2 « CL 3 / 2 v 

for some fixed C > 0. Note how the L 3 terms perfectly cancel. This finishes showing that 
||fi^l£) 2 ||i < e~ CL . As noted before, we may likewise develop a bound for ||i. 
Putting this together with (5.5) gives (5.3), which finishes the proof. □ 

Lemma 5.3. 

"Tilli 



in? Ih ► 0. 



L— >oo 

Proof. The proof of this result is the same as that of the previous lemma. Using the 
definition of Q 2 ^ and factorizing as in the above proof we get 

n 2 L = aq 2 l a* 

with 



„ -2mL+(A+A)L r 
&l(\~\) = t _ / dwdwe'^- -) 2 /8L+(«-+«.)L-4L3/3 

Jd 

■ Ai(-w + L 2 - X + m) Ai(-w + L 2 - X + m)l A ~ x<0 . 



Applying the change of variables w 1— > L 2 w and w 1— > L 2 w the kernel becomes 

\+~\)L/2-2mL r 

V8n JD 



r 7/2J\+\)L/2-2mL r 

n 2 L (X,X) = — — — — / dwdwe L f( - w '^Ai(L 2 (l-w) + m-X) 



where f(w,w) = ^-3 w ^ +(w + w) — i. Note the similarity with (5.4), the only difference 



•Ai(L 2 (l-u)) + m-A)l 

ote the simila 

being that in / we have a term — (w — w) 2 /8 instead of —(w + w) 2 /8. 

As in the above proof we need to bound ||0|||i, and to that end we split D into the same 
three regions D\, D2 and D' 2 . The operator restricted to D\ is easy to bound, exactly as 
before. On D 2 (and thus also on D' 2 ) we can repeat the same argument as before. The only 
diference is that, when we apply the change of variables w = 0+L~ 3 / 2 v and w = 4+L _3//2 {;, 
the function Hl{v) in the resulting integral in (5.7) is now multiplied by e 2vIj3/2 , coming 
from the difference between / and the function / defined after (5.4). This change does not 
affect the bound on the v integral (II and II in the above proof). It is straightforward to 
check that the rest of the proof is not affected either (note in fact that the only place where 
the definition of hi{v) is used is (5.10), and the approximation there is still valid). □ 

References 

[ACQ11] G. Amir, I. Corwin, and J. Quastel. Probability distribution of the free energy 

of the continuum directed random polymer in 1 + 1 dimensions. Comm. Pure 

Appl. Math. 64.4 (2011), pp. 466-537. 
[AGZ10] G. W. Anderson, A. Guionnet, and O. Zeitouni. An introduction to random 

matrices. Vol. 118. Cambridge Studies in Advanced Mathematics. Cambridge: 

Cambridge University Press, 2010, pp. xiv+492. 



CONTINUUM STATISTICS OF THE AIRY 2 PROCESS 



11 



[AS64] M. Abramowitz and I. A. Stegun. Handbook of mathematical functions with for- 
mulas, graphs, and mathematical tables. Vol. 55. National Bureau of Standards 
Applied Mathematics Series, 1964, pp. xiv+1046. 

[BFP08] F. Bornemann, P. Ferrari, and M. Prahofer. The Airyi Process is not the Limit 
of the Largest Eigenvalue in GOE Matrix Diffusion. J. Stat. Phys. 133 (3 2008), 
pp. 405-415. 

[BFS08] A. Borodin, P. L. Ferrari, and T. Sasamoto. Large time asymptotics of growth 
models on space- like paths. II. PNG and parallel TASEP. Comm. Math. Phys. 

283.2 (2008), pp. 417-449. 

[BR01] J. Baik and E. M. Rains. Symmetrized random permutations. In: Random matrix 
models and their applications. Vol. 40. Math. Sci. Res. Inst. Publ. Cambridge: 
Cambridge Univ. Press, 2001, pp. 1-19. 

[BS02] A. N. Borodin and P. Salminen. Handbook of Brownian motion — facts and for- 
mulae. Second Edition. Probability and its Applications. Birkhauser Verlag, 
2002, pp. xvi+672. 

[CH11] I. Corwin and A. Hammond. Brownian Gibbs property for Airy line ensembles. 

2011. arXiv:1108.2291. 
[Dys62] F. J. Dyson. A Brownian-motion model for the eigenvalues of a random matrix. 

J. Mathematical Phys. 3 (1962), pp. 1191-1198. 
[Erd56] A. Erdelyi. Asymptotic expansions. New York: Dover Publications Inc., 1956, 

pp. vi+108. 

[Fei09] T. Feierl. The Height and Range of Watermelons without Wall. In: Combinatorial 
Algorithms. Vol. 5874. Lecture Notes in Computer Science. Springer Berlin / 
Heidelberg, 2009, pp. 242-253. 

[Fis84] M. E. Fisher. Walks, walls, wetting, and melting. J. Stat. Phys. 34 (5 1984), 
pp. 667-729. 

[FMS11] P. J. Forrester, S. N. Majumdar, and G. Schehr. Non-intersecting Brownian 
walkers and Yang-Mills theory on the sphere. Nucl. Phys. B 844.3 (2011), pp. 500 
-526. 

[FNH99] P. J. Forrester, T. Nagao, and G. Honner. Correlations for the orthogonal-unitary 
and symplectic-unitary transitions at the hard and soft edges. Nucl. Phys. B 

553.3 (1999), pp. 601 -643. 

[FS05] P. L. Ferrari and H. Spohn. A determinantal formula for the GOE Tracy- Widom 

distribution. J. Phys. A 38.33 (2005), pp. L557-L561. 
[Joh03] K. Johansson. Discrete polynuclear growth and determinantal processes. Comm. 

Math. Phys. 242.1-2 (2003), pp. 277-329. 
[Mac94] A. M. S. Macedo. Universal Parametric Correlations at the Soft Edge of the 

Spectrum of Random Matrix Ensembles. Europhysics Letters 26.9 (1994), p. 641. 
[MQR11] G. Moreno Flores, J. Quastel, and D. Remenik. Endpoint distribution of directed 

polymers in 1 + 1 dimensions. 2011. arXiv: 1 106 . 2716. 
[PS02] M. Prahofer and H. Spohn. Scale invariance of the PNG droplet and the Airy 

process. J. Stat. Phys. 108.5-6 (2002), pp. 1071-1106. 
[PS11] S. Prolhac and H. Spohn. The one-dimensional KPZ equation and the Airy 

process. J. Stat. Mech. Theor. Exp. 2011.03 (2011), P03020. 
[QR12] J. Quastel and D. Remenik. Regularity and continuum statistics of the Airy\ 

process. 2012. arXiv: 1201 .4709. 
[RS10] J. Rambeau and G. Schehr. Extremal statistics of curved growing interfaces in 

1+1 dimensions. EPL (Europhysics Letters) 91.6 (2010). 
[RS11] J. Rambeau and G. Schehr. Distribution of the time at which N vicious walkers 

reach their maximal height. 2011. arXiv:1102 . 1640. 
[Sas05] T. Sasamoto. Spatial correlations of the ID KPZ surface on a flat substrate. 

Journal of Physics A: Mathematical and General 38.33 (2005), p. L549. 



CONTINUUM STATISTICS OF THE AIRY 2 PROCESS 



15 



[Sim05] B. Simon. Trace ideals and their applications. Second. Vol. 120. Mathematical 
Surveys and Monographs. American Mathematical Society, 2005, pp. viii+150. 

[TW96] C. A. Tracy and H. Widom. On orthogonal and symplectic matrix ensembles. 
Comm. Math. Phys. 177.3 (1996), pp. 727-754. 

[VS10] O. Vallee and M. Soares. Airy functions and applications to physics. Second. 
London: Imperial College Press, 2010, pp. x+202. 

(I. Corwin) Courant Institute, 251 Mercer St., Room 604, New York, NY 10012, USA 
E-mail address: corwin@cims.nyu.edu 

(J. Quastel) Department of Mathematics, University of Toronto, 40 St. George Street, 
Toronto, Ontario, Canada M5S 2E4 

E-mail address: quastel@math.toronto.edu 

(D. Remenik) Department of Mathematics, University of Toronto, 40 St. George Street, 
Toronto, Ontario, Canada M5S 2E4 

and Departamento de Ingenieri'a Matematica, Universidad de Chile, Av. Blanco Encala- 
da 2120, Santiago, Chile 

E-mail address: dremenik@math.toronto.edu 



